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Abstract 

The dependence of the pseudoscalar meson mass and decay constant is compared 
to one-loop Partially Quenched Chiral Perturbation Theory (PQChPT) in a numer- 
ical simulation with two light dynamical quarks. The characteristic behaviour with 
chiral logarithms is observed. The values of the fitted PQChPT-parameters are in 
a range close to the expectation in continuum in spite of the fact that the lattice 
spacing is still large, namely a ~ 0.28 fm. 



1 Introduction 

In numerical Monte Carlo simulations of QCD Chiral Perturbation Theory (ChPT) 
P is often used to guide the extrapolation to the physical values of the three light 
quark masses (m u ~ and m s ). In this procedure not only the lattice gauge theory 
results are established but also useful information is obtained about the values of 
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the Gasser-Leutwyler parameters of ChPT. In fact, recently several groups explored 
this possibility in quenched and unquenched simulations with Wilson-type [3] 
and staggered [1J|S] quarks. (For a review see 

In order to achieve small systematic errors the simulations themselves have to 
be performed in a range of quark masses where the applied one-loop (NLO) ChPT- 
formulas give a good approximation. In particular, the characteristic chiral loga- 
rithms have to appear in the quark mass dependence of different physical quantities. 
This applies both to original ChPT as well as to PQChPT [3 IE]- However, in most 
recent simulations - especially with Wilson-type quarks - this condition is not ful- 
filled because they are performed in the range m u ^ > \m s . Estimates based on 
present knowledge of the ChPT parameters indicate (see, for instance, that at 
least m u> d < jrn s — gm s has to be reached. (See also the summary of the panel 
discussion at the Boston Lattice Conference |1D|.) Another requirement is that the 
virtual effects of the s-quark also has to be taken into account by simulating with 
three light dynamical quarks. 

Since the dynamical quarks in most unquenched simulations do not satisfy the 
above bound, it is not surprising that the chiral logarithms have not been observed 
[TT1 IP21 IP31 ITU iTol ITS] . This was the main motivation of our collaboration to 
start exploring the possibility of simulating QCD with light quarks in the range 
m ud <im s |ni[TBlliniEniEI]- In these simulations we use the two-step multi- 
boson (TSMB) algorithm for dynamical fermions [22] and consider for the moment 
N s = 2 dynamical "sea" quarks. The case of N s = 3 is also under study [23] • Our 
first simulations were oriented towards the investigation of simulation costs as a 
function of the quark mass and were performed on modest size lattices (typically 
8 3 ■ 16) with lattice spacings of the order a ~ 0.27 fm - a value where continuum 
behaviour is not necessarily expected. Therefore, it came to us as a surprise that 
plotting the pseudoscalar ("pion") mass and decay constant as a function of the 
quark mass (in the form suggested by |24] 125]) the chiral logarithmic behaviour has 
been qualitatively displayed [El E01 • 

Encouraged by this result we picked out a point with m u d — \m s and performed 
a high statistics run on 16 4 lattice in order to study the dependence on the valence 
quark mass in a sufficiently large physical volume. This has the advantage that by 
taking ratios of the masses and decay constants the Z-factors of renormalization 
cancel. This removes an uncertainty in |19l I2U] where the Z-factors have been 
neglected by setting them to Z = 1. In our analysis of simulation data we applied 
PQChPT for Wilson lattice fermions [201 which take into account leading lattice 
artefacts of (9(a). 

The plan of this paper is as follows: in the next section the one-loop PQChPT 
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formulas for Wilson lattice fermions will be recapitulated. In section 3 the numerical 
simulation data will be analyzed and discussed. 



2 PQChPT formulas 

Our analysis of the valence quark dependence of the pseudoscalar mass (m n ) and 
decay constant (/„■) is based on the one-loop PQChPT formulas for the Wilson 
lattice action as derived in [201- Instead of the quantities with dimension mass- 
square xa and pa of ref. jSHj we prefer to use the dimensionless ones 



2B m q 2W acsw n \ 

XA = 72— > PA = 72 • C 1 ) 

Jo Jo 
Here m q is the quark mass, a the lattice spacing, Bq and Wq are parameters of 
dimension mass and (mass) 3 , respectively, which appear in the leading order (LO) 
chiral effective Lagrangian, csw is the coefficient of the 0{a) chiral symmetry break- 
ing term and /o is the value of the pion decay constant at zero quark mass. (Its 
normalization is such that the physical value is /o — 93MeV.) In ref. the case 
of three non-degenerate quark flavours is considered. Here we consider a general 
number N s of equal mass sea quarks. 

The next to leading order (NLO) PQChPT formula for the pion decay constant 
is in this case: 

=1 ~ 12871-2 i ^ XA + Xs + PA + Ps ^ log Q( Xj4 + XS + PA + Ps ^j 

+ {XB +XS+PB + PS) log Q(XB + XS + PB+ Ps)^j } 

+ 64iV 7T 2 i XA + XB + PA + PB ~ 2xs ~ 2ps + ( XB -XA + Pb - Pa)' 1 

■ 2(XA + Pa){xb + Pb) - (xs + Ps)(XA + XB + PA + Pb) log ( XA - PA ) I 

\XB + Pb) J J 

+2l 5 ( X A + Xb) + 2W 5 (p A + p B ) + 4iV s Z 4 X5 + ^N s W iPs . (2) 

Here A and B denote generic quark indices: S will be the label for the sea quarks 
V for valence quarks. For the pion mass squared we have: 

m AB l ( , , n . 1 (XA + XB + PA + Pb) 

- -AXA + XB + PA + Pb) + 
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{Xb + Pb){xs -Xb + Ps - Pb) log(xB + Pa) } 



{ 



+AN s {2i 6 ~ U)xs{xa + Xb) + 2 (2L 8 - L 5 ){ X A + Xb? 
+4N S (W 6 - L A ) Xs {pA + Pb) + 4N S (W 6 - W 4 )ps(XA + Xb) 
+2(2W 8 -W 5 -L 5 )( XA + Xb)(pa + Pb) ■ (3) 
The NLO parameters L k and W k are related to L k and W k in ref. [26 by 

L k = L k - c k log(/ 2 ) , W k = W k -d k log(/ 2 ) , (4) 

where the coefficients of the logarithms are given by 

1 N 8 (iV 2 + 2) (ATf-4) 

4 256vr 2 ' 5 256^ 2 ' 6 5127V 2 tt 2 ' 8 512iV s 7r 2 ' 1 > 

respectively, 

j - -J— - - ^ + 2 ) w - (^ 2 ~ 4 ) ™ 

4 " 256vr 2 ' 5 " 256tt 2 ' 6 " 256iV|7r 2 ' 8 " 256iV s 7r 2 " W 

The relations in @ have the unpleasant feature that logarithms of a dimensionful 
quantity appear. One can avoid this by introducing 

L' k = L k - c k log(^ 2 ) , W' k = W k - d k log(/i 2 ) , (7) 

where [i is the mass scale introduced by dimensional regularization. Since L k and 
W k depend on [i the choice of it in the logarithm is natural. In terms of L' k and WL 
we have 

L k = L' k -c k log(^pj , w k = W^-d k \og(Jpj . (8) 

Note that the NLO parameters a k in ref. are related to L' k by 

a k = 128ir 2 L' k . (9) 

The universal low energy scales A% t 4 in ref. |241 125j can be expressed, in the case of 
N s = 2, by the following combinations of the coefficients L k : 

^log| = 2l 4 + I 5 . (10) 

In this paper we are interested in the valence quark mass dependence of /„■ and 
m 2 for fixed sea quark mass parameter \s- Therefore it is natural to introduce the 
ratios of the other mass parameters to \s'- 

t - Xv - Ps A _ PV PV 

Xs Xs Ps VXS 



Once relations (|TT|) are substituted in (JSJ-Q, the logarithmic dependence on xs 
can be absorbed in the NLO parameters if one introduces 

L S k = L k -c k log(xs) = L' k - c fe log (j^XS^J , 

W Sk = W k -d k \og( XS ) = W' k -d k \og{^ X ^j ■ (12) 

Let us note that the argument of the last logarithms here can also be written as 

/ 2 2B m qS , , 

— XS = 2 • ( 13 ) 

In this paper we keep the sea quark mass (xs) fixed and vary the valence quark 
mass (xv = £,Xs)- Expanding the ratio of decay constants up to first order in the 
one-loop corrections one obtains 

Rf vv = lpL = i + 4(^- l)xsL S5 

JSS 

NsXS , t , o m 1 + £ + 2 V , N sXS n . \ I / 1 i \ f- [A \ 

~ 64^ + ^ + 2r ?) lo S 2 32^ n ' g ^ ^ ' ^ ' 

and similarly 

Rhs ^ |f = 1 + 2({ - l) XS L Si + j^-jK - 1) - ^(1 + ,) log i±2 

-^(1 + « + 2,) log + + „) log(1 + „) . (15) 

In case of the ratios of m 2 we expand up to first order in the "correction" which 
now also includes the 0{a) terms of the tree- level expressions: 
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+8£(£ ~ Vxs(2L S8 - L S5 ) + 8N S (£ - l)rixs(Ls4 - W S6 ) 
XS r (e-l)(^ + ? ?)--^ T e(l + 2r ? )log(l+r / ) 



and 
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+ TFrh2^ 2 ~ £ ~ * + 3r >0 iog(e + v) , (i6) 
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Rm vs = ^ = ^(1 + £ + rj - rjO 
m SS 2 

+2(£ + 1)(£ - l)xsC2L S8 - L S5 ) + 4iV s (e - l)vXs{L S 4 - W se ) 
Xs 



32iV s vr 2 



(C + 1)(1 + 27?) log(l + r/) 



+-^(^ 2 + £ + V + 3r/0 log(£ + r,) . (17) 



In these expressions it is assumed that the 0{a) mass terms p are the same for 
valence quarks and the sea quark, namely py = ps- This is the case if only the 
hopping parameter k is changed. Changing py = CvXs can be investigated by 
changing the Wilson-parameter r in the Wilson fermion action, too. 

Up to now we considered the valence quark mass dependence for unchanged sea 
quark masses. Let us remark that the formulas for the sea quark mass dependence 
can also be written in a similar form as eqs. (|14 j) -(|17 j) . In this case it is advantageous 
to fix a reference sea quark mass XR ( see 0) and introduce the variables 

_ XS _ PR _ PS PS ,, Q , 
a = — , ijj = — , t = — = . (18) 

XR XR PR UXR 

Instead of the NLO parameters in Q12j) the appropriate ones are then obviously 

L Rk =L' k -c k log(^ X R^ , W Rk = W' k -d k \og{^ X ^j • (19) 

3 Numerical results 

We performed Monte Carlo simulations with N s = 2 degenerate sea quarks on a 16 4 
lattice at (3 = 4.68, k = 0.195 and investigated the valence quark mass dependence 
at k = 0.1955, 0.1945, 0.1940, 0.1930, 0.1920. The statistics corresponds to 1180 
gauge field configurations. The error analysis was based on the linearization method 
[27] . Since r^/a = 1.76(6) the lattice spacing is a ~ 0.28 fm. This means that the 
physical lattice extension is rather large: L ~ 4.5 fm. The value of the quark mass 
parameter is given by am w = 0.519(1) as M r = (rQm w ) 2 ~ 0.83. This is about \ of 
the value of M r for the strange quark M^ trange ~ 3.1. 

The ratios in eqs. (jl4[) -(|17 |) as a function of the quark mass ratio (£) depend 
on five parameters, namely xs, Vi ^S5, (2-^58 — -^ss) and (L$a — Wse)- With our 
choice of the valence hopping parameters and with our statistics most of the multi- 
parameter fits were unstable therefore our analysis is based one a sequence of single 
or double parameter fits. The stability of the multi-parameter fits can be improved 
by optimizing the choice of valence quark mass values, which we did not exploit this 
time. 

A very useful quantity is the double ratio of decay constants ^1] which does not 
depend on any of the NLO coefficients. In other words there one can see the chiral 
logarithms alone. The NLO formula is: 

W = 7 ^ = l + ^- f (C-l)-^ T (l + r/)logf^ . (20) 
Jvvfss 32N s tt z 32N s tt z l + r) 

Because of the 0(a) contributions this has two parameters: xs an d r] = ps/xs- F° r 
performing two-parameter fits some timeslice distance pairs were chosen and the fits 
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Valence quark mass dependence of RRf 
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Figure 1: The valence quark mass dependence of the double ratio of 
pion decay constants RRf. Besides the "fit" the two other curves 
show the 0(a) contribution ("eta") and the physical contribution 
obtained at rj — ("chi"). 



of the pion- and quark mass were taken from them. Useful choices are, for instance, 
4-5, 4-6 or 5-6. (The way physical quantities were obtained has been described in 
detail in [T7] ) . The result of the two-parameter fit was (see figure [""J : 

Xs = 9.2±2.5, n = 0.14 ±0.30 . (21) 

The value of xs has the right order of magnitude. Indeed, from the axial Ward iden- 
tity we obtain rom q s = 0.06 Z" 1 . Using this and the phenomenological estimates 
[2*5] ro/o = 0.23, rn-Bn = 7.0, where the value of Bq refers to the MS scheme at 
[i = 1 GeV, we deduce ^ stmate ~ 16 Z" 1 . Here Z q is an unknown Z- factor relating 
the bare lattice quark mass to the renormalized one at 2 GeV, which is typically 
of 0(1). Another estimate can be obtained by using the tree-level ChPT formula 

x estimate _ M r /(r f ) 2 ~ 15.7. 

As one can see from (|21|). the value of the ratio n = ps/xs has a large error in 
the fit. The reason is that the functional form of the terms due to the leading 0{a) 
lattice artefacts is not very different from those coming from the non-zero quark 
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Valence quark mass dependence of LRm 
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Figure 2: The valence quark mass dependence of the linear combi- 
nation of pion mass-squared ratios LRm. Besides the "fit" the two 
other curves show the O(a) contribution ("eta") and the physical 
contribution obtained at n = ("chi"). 



mass. Nevertheless, they are different and an optimized choice of ^-values plus high 
statistics allowing for stable multi-parameter fits makes the distinction of the two 
sorts of chiral symmetry breaking possible. In the rest of this paper we only consider 
one- and two-parameter fits with fixed r] in the range given by l|21|). 

After determining xs and n = ps/xs from the double ratio RRf one can fit the 
other ratios to obtain estimates of the NLO coefficients. A nice linear combination 
of mass-squared ratios is: 

LRm = 2Rm vs - Rm VV = 1 - 4(£ - l) 2 Xs(2L S8 - L S5 ) - -^^(C - 1)(£ + rj) 

- tM^ {1 + 2l]) log(1 +v) ~ " 2i " 2l]) log(e + v) ■ (22) 

This has only a single new parameter (2L$8 — Lgs) and the statistical errors are 
small, therefore one can also perform a two-parameter fit of xs an d (2Lss — Lss) 
with the result (n = 0.05 fixed, errors in last digits given in parentheses): 

2L 58 - L S5 = -0.00203(5) , 2a 8 - a 5 = 0.85(6) , X s = 5.2(1.1) . (23) 
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Valence quark mass dependence of RfVV 
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Figure 3: The valence quark mass dependence of the ratio of pion 
decay constants Rfvv- Besides the "fit" the two other curves show 
the O(a) contribution ("eta") and the physical contribution obtained 
atr) = Q ("chi"). 



oik denote the NLO parameters in Q taken at the renormalization scale /x = 47r/o- 
Fixing both \s = 10.0 and rj = 0.10 gives: 

2L S8 ~ L S5 = -0.00177(3) , 2a 8 - a B = 0.58(3) . (24) 

As figure [21 shows, with these parameters the last point is not perfectly fitted. A 
perfect fit is obtained with the parameters in (|23|) . 

The value of L55 can be determined from Rfvv- {Rfvs gives very similar re- 
sults.) In this case the statistical errors are larger, therefore only a single parameter 
fit is useful. The result for fixed \s = 10.0 and t/ = 0.10 is (see figure 

L 55 = 0.0034(1) , a 5 = 1.55(24) . (25) 

The errors given in (j23|) -(J25 J) are the ones for the specified values of fixed parameters. 
The overall error is, of course, larger - as one can see, for instance, by comparing 
(|23j) and (|24*|) . The values of (2a§ — "5) and 0:5 are somewhat larger than the results 
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of UKQCD 0: (2a 8 -a 5 ) = 0.36±0.10 and a 5 = 1.22±0.11 (only statistical errors 
quoted). 

The double ratio of the pion mass squares |13| 



can be used to determine the fifth parameter (L54 — Wse)- In this case a single 
parameter fit with fixed xs = 10.0 and r\ = 0.10 gives (L54 — Wse) = 0.00358(6). 

The conclusion of this paper is that - once the quark masses are small enough - 
the qualitative behaviour of the low energy chiral effective theory with chiral loga- 
rithms is present even on coarse lattices. Since here ratios of pion masses and decay 
constants are considered the Z-factors of renormalization cancel, therefore the un- 
certainty about their /3-dependence, which can influence the results of ^UEI]> is 
removed. The coefficients of the observed chiral logarithms and the fitted values of 
the Gasser-Leutwyler coefficients are close to expectation. This qualitative agree- 
ment of the results of a numerical simulation with (PQ)ChPT is quite satisfactory 
but for a quantitative determination of NLO ChPT parameters one has to perform 
extrapolations to a — > and m q — > 0. Since 0{a) effects are taken into account in 
the analysis by the Rupak-Shoresh effective Lagrangian, the continuum limit will 
be reached asymptotically at the rate 0(a 2 ). 

The computations were performed on the APEmille systems installed at NIC 
Zeuthen, the Cray T3E systems at NIC Jiilich and the PC clusters at DESY Ham- 
burg. We thank Hartmut Wittig for helpful discussions on various topics and for 
careful reading of the manuscript. I.M. profited from enlightening discussions on 
the Asia-Pacific Mini- Workshop on Lattice QCD especially with Shoji Hashimoto, 
Yoichi Iwasaki, Yusuke Namekawa and Steve Sharpe. 
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